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ABSTRACT 

A radiative transfer algorithm in the solar wavelengths 
for the NOAA POES AVHRR and GOES Imager is proposed for the 
cloud-free, marine atmosphere. The algorithm combines 
linearized, single-scattering theory with an estimate of bi- 
directional surface reflectance. Phase functions are 
parameterized using an aerosol distribution model and the 
ratio of radiance values measured in channels 1 and 2 of the 
AVHRR. Retrieved satellite aerosol optical depth is 
compared to airborne sunphotometer data and values derived 
from aerosol particle size distributions collected during 
the Tropospheric Aerosol Radiative Forcing Observational 
Experiment (TARFOX) in July 1996. Error in the satellite 
derived values from the AVHRR originates in error in 
modeling aerosol size distributions, corresponding phase 
function parameterization and treatment of specular surface 
reflectance. Extension of the algorithm to the GOES Imager 
provided results consistent with the AVHRR. 
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I. 



INTRODUCTION 



Atmospheric aerosols play a major role in the Earth's 
radiation budget. Aerosols contribute to both a direct and 
indirect change in the Earth's albedo which may lead to 
changes in the heating of the Earth and its atmosphere. 
Variations in natural and anthropogenic aerosols directly 
alter the scattering of incoming shortwave radiation to 
space. Particles that become cloud condensation nuclei 
create variations in local cloud characteristics which 
indirectly affect solar insolation. Charlson et al . (1992) 
report that shortwave radiative forcing of the Earth's 
climate due to an estimated increase in anthropogenic 
sulfate optical depth of 0.04 is enough to offset longwave 
forcing by greenhouse gases. Recent research efforts have 
focused on anthropogenic sources of aerosols in coastal 
regions and their radiative effect on the Earth's climate. 

For the U.S. Navy, characterization of aerosol 
radiative properties in the coastal zone is important to the 
design, planning, and operation of electro-optical weapons 
and sensor systems near coastal boundaries. Aerosols 
scatter and degrade electro-optical system performance by 
scattering visible and near-infrared energy used by these 
systems . 

In situ measurements of both aerosol radiative and size 
distribution characteristics using instruments such as 
sunphotometers, spectrometers, and radiometers, provide 
highly accurate, but spatially and temporally limited 
measurements. Satellite radiometers offer an alternative 
method to measure aerosol radiative effects that cover large 
areas . 
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To be useful for this purpose, a satellite radiometer 
must possess high spatial and radiometric resolution in its 
visible and/or near-infrared channels. Until recently, the 
only operational U.S. satellite system with fine enough 
spatial and digital resolution to indirectly measure aerosol 
radiative properties was the NOAA Advanced Very High 
Resolution Radiometer (AVHRR) carried aboard the NOAA Polar- 
orbiting Operational Environmental Satellites (POES) . With 
the POES series satellites, retrievals are limited to one or 
two passes per day over a given mid latitude area due to the 
nature of the polar orbit. Starting in late 1994, NOAA 
launched a series of two new geostationary weather 
satellites, the Geostationary Operational Environmental 
Satellites or GOES. GOES included an improved visible 
spectrum radiometer (GOES Imager) with similar radiometric 
resolution to the AVHRR. The primary advantage of the GOES 
series radiometer is its capability of imaging an area 
approximately every 15 minutes. This enhanced temporal 
resolution potentially offers insight into short-term, 
large-scale aerosol variations previously not available. 

Most of the satellite aerosol retrieval techniques to 
date have focused on the AVHRR. Unique to the AVHRR among 
operational meteorological satellites is a two solar channel 
capability that allows better characterization of phase 
scattering effects; this capability is not available using 
the single solar channel GOES Imager. Earlier work by 
Durkee et al . (1991) and Rouault and Durkee (1992) exploits 
this capability of the AVHRR and provides the basic 
framework for the aerosol retrieval technique used in this 
study. 

Until recently, the ability to validate optical depth 
retrieval results has been limited by the lack of field data 
to provide closure on comparable measurements. In July 
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1996, the Tropospheric Aerosol Radiative Forcing 
Observational Experiment (TARFOX) was conducted to attempt 
to provide this closure. 

The objectives of this thesis are threefold: 

- Develop an improved linearized, single-scattering 
atmospheric aerosol optical depth retrieval algorithm for 
use by both NOAA POES AVHRR and GOES Imager systems based on 
the approach of Durkee et al . (1991). 

Validate aerosol optical depth retrievals using 
TARFOX field observations of optical depth and aerosol size 
distributions . 

- Demonstrate temporal utility of GOES aerosol optical 
depth retrievals. 

Chapter II describes the basic radiative transfer 
theory used in the satellite optical depth retrievals. 
Chapter III outlines the TARFOX data set and basic 
instrumentation used. Chapter IV describes aerosol optical 
depth retrieval procedures. Chapter IV provides results. 
Chapter VI discusses final conclusions and recommendations. 
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II. RADIATIVE TRANSFER THEORY 



Radiative transfer theory provides the basis for 
methods used to characterize aerosol properties from 
satellite remote measurement (see Liou, 1980, for a complete 
treatment of the theory) . In a cloud-free, marine 
environment, the shortwave, solar radiation measured by a 
satellite radiometer is primarily the result of scattering 
by both the molecular constituents of the atmosphere 
(Rayleigh scattering) and larger suspended aerosol (Mie 
scattering) . Absorption plays a small role in the 
attenuation of incoming solar radiation in the visible (0.4- 

0.7 |im) and near-infrared (NIR) spectrums (0.7-0. 9 Jim). 
Corrections to solar irradiance of less than 5% for ozone 
absorption in the upper atmosphere are applied; aerosols are 
assumed to be non-absorbing. At some solar wavelengths, 
corrections for water vapor absorption on the order of 2-3% 
of measured radiance are also applied. In the absence of 
Sun glint, reflectance from the ocean surface is also small. 
Contributions to satellite measured radiance due to sea 
surface foam and subsurface reflectance can be accounted for 
with empirical measurements. Specular reflectance is 
estimated using Fresnel coefficients. Figure 2.1 
illustrates this radiative transfer process. This chapter 
describes the basic theory used in the methods employed by 
this study to obtain aerosol optical properties. All 
further discussion pertains to radiative transfer in a 
cloud-free, marine environment. 

A. EXTINCTION, SCATTERING, AND OPTICAL DEPTH 

Radiative extinction includes both scattering and 
absorption by atmospheric constituents; extinction is highly 
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wavelength dependent. Equation 2.1 defines the extinction 
coefficient: 



o ext = J 7rr 2 Q ext (m, r)n(r)dr (2.1) 

0 

where r is particle radius, nr 2 is particle cross-sectional 
area, Q ext (m,r) is the extinction efficiency factor, m is 
the complex index of refraction, and n(r)is the number of 
particles for a given radius. Q ex t( m / r ) is a function of 
both composition and size of a particle and describes the 
effects of both scattering and absorption due to the 
interaction of a particle with radiative energy of a 
specified wavelength. Changes in the size, composition or 
distribution of constituents or suspended particles in the 
atmosphere directly affect the amount of extinction 
observed. Therefore, measurement of extinction can 

inversely lead to knowledge of the characteristics of the 
atmosphere's particulate distribution. 

Scattering of solar radiation dominates other radiative 
transfer factors in the visible and NIR. Due to the near 
uniformity of molecular constituents both spatially and 
temporally for a given region, Rayleigh (molecular) 
scattering can be adequately calculated (Durkee et al . , 
1991) . Removal of Rayleigh scattered radiance effects from 
satellite retrieval methods allows for quantification of 
scattering effects by suspended aerosols only. The effect 
of scattering due to aerosols is approximated using Mie 
theory for spherical particles. For the cloud-free marine 
environment, the extinction coefficient is well represented 
by the scattering coefficient: 
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scat 



a 



= / 7ir 2 Q scat (m, r)n(r)dr 



( 2 . 2 ) 



where Q scat (m,r) is the scattering efficiency factor, 

representing the ratio of total energy scattered in all 

directions to incident energy. 

By integrating the extinction coefficient in the 

vertical through the atmosphere, optical depth is obtained: 

H H 

6 = J tfext dz * 1 a scat dz (2.3) 

0 0 

Because satellite radiometers measure the radiative 
properties of the entire atmospheric column, quantification 
of total column optical depth is the objective of the 

satellite retrieval technique used in this study. 

B. RADIATIVE TRANSFER SOLUTION 

Equation 2.4 provides a general form for radiative 
transfer for a given solar wavelength in a plane parallel 
atmosphere is defined as (after Liou (1980)): 






dL t (0; |i, <|>) 
dS 



©or 

L t (5, Cl) - — ! L t (5, C1)&C1, Cl )dQ 
4n 4n 

m - 5/ 

- — 7cF 0 p(f2 - Cl Q )e /Ho 



(2.4) 



where n is the cosine of satellite zenith angle (9) , [i is 

the cosine of Sun zenith angle (0o) , (p is the relative 
azimuth between the satellite and Sun, L t is diffuse 
radiance (W/m 2 )im sr) , £> is the solid angle defined by 0 and 

cp, 0) o is the single scattering albedo, p is the scattering 
phase function, and F 0 is the incoming solar radiance at the 
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top of the atmosphere. The single scatter albedo (co 0 ) is 
defined as the ratio of the scattering coefficient to the 
extinction coefficient. The terms on the right hand side of 
Equation 2.4 represent the total radiative energy in a beam, 
the energy scattered into the beam due to multiple 
scattering, and the energy scattered into the beam due to 
single scattering, respectively. 

For atmospheres with small optical depths such as the 
cloudless, marine atmosphere, contribution to Equation 2.4 
by multiple scattering is negligible. Solution to Equation 
2.4 without multiple scattering is known as the single 
scattering approximation (after Liou (1980)): 



Lt 



<°oHo F o 

ip + Ho) 





(2.5) 



where \j/ s is the scattering angle. 

After accounting for ozone, Rayleigh and ocean surface 
effects, Durkee et al . (1991) show that in this environment 

8 a is small enough to reduce Equation 2.5 to: 

L a = ^^p(vFs) 5 a (2.6) 

where the subscript “a” denotes aerosol related quantities. 

Solution of Equation 2.6 for aerosol optical depth (8 a ) 
requires knowledge of the measured radiance at the satellite 
radiometer (L a ) due to aerosol scattering at a given 
wavelength, the single scatter albedo, the incoming solar 
radiance (F 0 ) , the satellite zenith angle (represented by 

(a) , and the scattering phase function (P) . Since typical 
marine aerosol are weak absorbers, the scattering and 
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extinction coefficients are approximately equal for this 
environment, and the single scatter albedo is approximately 



satellite-Earth geometry. The other terms are described in 
detail below. 

1. Aerosol Radiance (L a ) 

Equation 2.7 describes the aerosol scattered radiance 
measured by a satellite radiometer through a simple, linear 
expression (after Gordon and Clark, 1980) : 



where L r is radiance due to Rayleigh Scatter, L s is surface 
reflected radiance, L g is Sun glint radiance, and x is 
transmittance of the atmosphere. For satellite radiometer 
window channels in solar wavelengths, x can be approximated 
as one with less than 10% error. 

To remove ozone absorption effects in the upper 
atmosphere, a correction to L t is made by dividing L t by the 
slant-path transmittance due to atmospheric ozone as in 
Equation 2.8: 



one. 



The satellite zenith angle is determined by the 



I. 



' a 




(2.7) 



I, 




( 2 . 8 ) 



where 5 03 is ozone optical depth. 
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Radiance contribution associated with Rayleigh scatter 
is calculated using a two-stream model by Turner (1973): 



Lr = 



E 0 e 



- 8 , 



03 



(%„*%) 






K S r[p r K) + P'(V..)] + 

+ *•’ :;■(! - p) 5l 

[K [P r (v s + ) + P r (M/ s -)] - 



-8 r / 

H 1 - e /fl 



\^c 



1 + (l - p)5 r 



(5 r + n)e ^ - n 



(2.9) 



where E 0 is incoming solar irradiance, 5 r is Rayleigh optical 
depth, P r (M / s ) i s the Rayleigh scattering phase function 
(equal to 3/4 (l+cos2 (\j/ s ) ) where + indicates forward 

scattering and - indicates back scattering, and p is 
reflectivity of seawater. 

L s is not modeled. Ramsey (1968) , using Coastal Zone 
Color Scanner (CZCS) , demonstrated that surface reflectivity 
is less than 0.005 for red-visible wavelengths and zero for 
wavelengths greater than 0.7 pm, consistent with the 
approximation of the ocean as a blackbody. Surface 

reflectance due to foam on the sea surface and subsurface 
reflectance due to suspended plankton/chlorophyll are the 
major contributors to L s . Koepke (1984) provides a measure 
of foam reflectance in the visible spectrum based on wind 
speed from surface leaving radiance measurements ranging 
from near-zero in winds below 5 m/s to 0.02 in winds in 
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excess of 20 m/s. Ignatov et al . (1995) report subsurface 
reflectance at 0.63 Jim in the open ocean with chlorophyll 
concentrations less than 0.25 mg/m 3 on the order of 
0.001410.0006. For this thesis, the combined effect of foam 
and subsurface reflectance is estimated at 0.002 at 0.63 pm 

and 0.005 at 0.86 pm. Directional reflectance due to the 
variation in refractive index between the atmosphere and 
ocean is treated in the phase function term described below. 

Contamination by specular glint (L g ) of the surface of 
the ocean for low Sun angles can be a problem for certain 
satellite-Earth-Sun geometries. Special care must be taken 
to review a satellite image for Sun glint signatures and 
avoid optical depth retrievals in contaminated regions. L g 
is not modeled; in glint free areas L g is assumed to be 
zero. 



2. Scattering Phase Function (P) 

The scattering phase function describes the probability 
of energy scattering into a specific direction. The 
scattering phase function is dependent on the radiation 
wavelength and size, composition, and distribution of 
atmospheric constituents. For Rayleigh scattering, the 
phase function is well understood; this approximation is 
used in the Turner solution above (Equation 2.9). For 
larger particles, Mie theory is used to quantify the 
scattering phase function with wavelength, size distribution 
and refractive index as inputs. To accurately quantify the 
phase function, one must have knowledge of the aerosol size 
distribution and composition. This creates an ill-posed 
problem since it is precisely the aerosol characteristics 
desired as a result of this retrieval method. To 
parameterize the phase function, several methods may be 
used. A common approach incorporates a phase function 
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empirically fit using measured averaged size distributions 
and characteristics for a region. Another uses the 
sensitivity of the phase function to radiative differences 
between two (or more) wavelengths to parameterize the phase 
function. Chapter IV contains a detailed description of 
both methods. 

When solar radiation scatters due to molecular or 
aerosol interaction, the scatter can be described as both 
direct (directly towards the receiving sensor) or diffuse 
(scattering at other angles which may eventually scatter or 
reflect back to the sensor) . Based on Sun-Earth-satellite 

geometry, both a forward (0-90°) and back scattering angle 

(90-180°) are defined. Treatment of the direct scatter is 
accomplished by using the value of the phase function in the 
linear single scattering model corresponding to the back 
scattering angle. Without true multiple scattering in the 
model, effects of diffuse scatter reflecting off the ocean 
surface (bi-directional reflection) can be estimated using 
Fresnel reflection coefficients and the forward scattering 

angle (0-90°). An effective phase function is defined by 
Equation (2.10): 



Peff = P- + P+[ r n 0 + r n] (2.10) 

where P e££ is the effective phase function, P. is the value 
of the phase function at the back scattering angle, P. is 
the value of the phase function at the forward scattering 
angle, and r is the Fresnel reflection coefficients at both 

p 0 and p. Fresnel reflection coefficients are given in 
Equation (2.11): 
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r e 



0.5 



( 2 . 11 ) 



f sinte - e ^ 2 

vsin^ + d t )j 



^tan 

K tan(0 i + 0 t ) y 



where 0 t is the angle of incidence and 0 t is the angle of 

transmission. From Snell's law, 0 t = sin' 1 (sin (0 X ) /m) ) where 
m is the index of refraction of seawater (1.33). The 
effective phase function describes an estimate of the 
increase in probability of solar radiation received at the 
sensor due to two possible paths: (1) direct transmission of 

solar energy to the sea surface, air-ocean interface 
reflectance, and forward scatter off an aerosol, or (2) 
forward scatter off an aerosol, air-ocean interface 
reflectance, and direct transmission to the sensor. 



13 




Figure 2.1. Radiative Transfer in the cloudless, 
marine environment. On the left, a typical 
atmospheric temperature profile depicting the marine 
atmospheric boundary layer (MABL) is included. Water 
vapor and aerosols are assumed to be confined to the 
MABL. Various scattering paths of incoming solar 
radiation (F 0 ) from Sun to satellite are shown in the 
right panel. Radiance variations along these paths 
measured by satellite radiometers (L t ) are primarily 
caused by upper atmosphere ozone absorption, Rayleigh 
(L r ) and aerosol (L a ) scattering, and ocean surface 
reflectance (L s ) . 
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III. DATA 



In order to validate the optical depth retrieval method 
described in this study, reference data sets were chosen 
based on the availability of in situ measurements of aerosol 
distributions and/or optical depth measurements. To the 
greatest extent possible, in situ measurements are matched 
to satellite observations both spatially and temporally. 
This chapter will briefly describe the data sets and 
instrumentation used to collect the data used in this study. 

A. TARFOX DATA SETS 

The TARFOX field experiment was designed as a closure 
study to better understand the radiative forcing effects of 
aerosols. TARFOX was conducted in the continentally- 
influenced environment off the eastern coast of the United 
States near Wallops Island, Virginia from 10-31 July 1996. 
During TARFOX, a variety of aerosol conditions ranging from 
relatively clean to moderately polluted were observed and 
measured. In situ measurements were conducted by airborne 
platforms including the University of Washington's (UW) 
C-131A and the United Kingdom's (UK) Meteorological Research 
Flight C-130 (see Russell et al . , 1996). The UW C-131A was 
outfitted with NASA AMES Airborne Autotracking Sunphotometer 
for aerosol optical depth, a Passive Cavity Aerosol 
Spectrometer Probe (PCASP) for aerosol size distributions, 
and various meteorological observing instruments 
(temperature, dew point, pressure) . The UK C-130 carried a 
PCASP, a Particle Soot /Absorption Photometer (PSAP) to 
measure aerosol absorption, and various meteorological 
observing instruments. Satellite imagery collected included 
full resolution NOAA 14 AVHRR (5 channel) and GOES 8 Imager 
(channel 1 and channel 4) data. The TARFOX Operations 
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Summary (Whiting et al . , 1996) contains details of the field 
collection effort. 

B. INSTRUMENTS 

1. NOAA Advanced Very High Resolution Radiometer 
(AVHRR) 

The AVHRR instrument is a component of the NOAA Polar 
Orbiting Operational Environmental Satellite (POES) series 
satellites. Current operational POES include the NOAA 12 
and NOAA 14. These satellites are in Sun synchronous orbit 
(883 km) and provide two passes per day in the morning and 
evening, respectively. Due to current orbit design, NOAA 12 
does not provide adequate daylight morning imagery for use 
in the aerosol optical depth retrieval in the mid latitudes 
and was not used in this study. The AVHRR instrument 

measures radiant and solar-reflected energy from sampled 
areas of the Earth in five spectral bands with a sub- 
satellite resolution of 1.1 km. Table 3.1 lists the 

characteristics of the individual radiometer channels for 
the AVHHR . (Kidwell, 1995) 



Channel 


Band Widths (|im) 


1 (Visible) 


0.58 - 0.68 


2 (NIR) 


0.725 - 1.10 


3 (IR) 


3.55 - 3.93 


4 (IR) 


10.3 - 11.3 


5 (IR) 


11.5- 12.5 



Table 3.1. NOAA AVHRR Radiometric 
Channels 



Channels 1 and 2 are used in the optical depth 
retrieval. Channels 4 and 5 are used in manual cloud 
screening analysis. 
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All AVHRR channels are calibrated prior to launch. 
Channels 1 and 2 of the AVHRR have no onboard calibration 
systems. Post calibration methods for these channels have 
been developed by the NOAA/NESDIS Office of Research 
Applications based on the results of Rao and Chen (1995) . 
Calibration formulae are incorporated into the satellite 
image processing discussed in Chapter IV. 

2 . GOES Imager 

This instrument is a major component of the NOAA 
Geosynchronous Operational Environmental Satellite (GOES) 
series satellites GOES 8 and GOES 9. Operating in a 
geosynchronous orbit at 36,000 km over equatorial sub-points 
at 75 W (GOES 8) and 135 W (GOES 9) , this radiometer 
provides imagery from full-Earth disc images to small area 
scans due to a flexible scan system design. Like the AVHRR, 
the GOES Imager is a multiple channel radiometer which 
measures radiant and solar-reflected energy from the Earth. 
Unlike the AVHRR, spatial resolution is not the same for 
each channel. All GOES Imager channels are calibrated prior 
to launch. Similar to the AVHRR, the GOES Imager channel 1 
does not have an onboard calibration system. No post- 
calibration methods for this channel have been developed or 
are planned for development. Table 3.2 lists the spectral 
and resolution characteristics of the GOES Imager. Channel 
1 is used in the optical depth retrieval method described in 
this study. ( GOES-IJ/KLM SN03 Imager Data and Calibration 
Handbook , 1994, and GOES Calibration and Alignment Handbook 
for the Imager SN04 Instrument, 1994) 
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Channel 


Band Width (|im) 


Resolution (km) 


1 (Visible) 


0.55 - 0.75 


1 


2 (IR) 


3.80 - 4.00 


4 


3 (Water Vapor) 


6.50 - 7.00 


8 


4 (IR) 


10.20 - 11.20 


4 


5 (IR) 


11.50 - 12.50 


4 


Table 3.2. GOES Imager Radiometric Channel 


Ls 



Only imagery from GOES 8 was used in this study due to 
the geographic location of TARFOX. 



3. NASA AMES Airborne Autotracking Sunphotometer 
By tracking the Sun, the NASA AMES Airborne 

Autotracking Sunphotometer measures the relative intensity 
of the direct incoming solar radiation in multiple spectral 
channels. This information is converted to optical depth. 
By flying an aircraft at low altitudes in cloud-free regions 
within the atmospheric boundary layer, total column optical 
depth can be well approximated. Table 3.3 lists the 

spectral bands for the six channels of the instrument. 
Instrument calibration is performed by NASA AMES prior to 
flight; the instrument is designed to maintain calibration 
within 1% during operation. Resolution of optical depth is 
on the order of 0.01. (Matsumoto et al . , 1987) 
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Channel 


Spectral Band (urn) 


1 


0.380 


2 


0.451 


3 


0.525 


4 


0.860 


5 


0.945 


6 


1.021 



Table 3.3. NASA AMES Airborne 
Autotracking Sunphotometer Spectral 
Characteristics 



During TARFOX only channels 1-3 and 6 were reported. 

4. Passive Cavity Aerosol Spectrometer Probe (PCASP) 
The PCASP is an airborne instrument designed to size 

aerosol distributions in the 0.1 (im to 3.0 |im range. The 
PCASP employs a He-Ne (632.8 nm) laser to sample aerosol 
size distribution from an aerodynamically focused jet which 

restricts the particle flow to a 150 jj.m diameter stream. A 
combination ref lecting-refracting optical system collects 
light scattered from the laser beam by the aerosol in the 
jet and converts the resultant signal into particle size 
distribution. Distribution is segmented into 15 channels 

with a minimum size resolution of 0.02 Jim in the smallest 
size channels, progressively weighted towards a maximum of 
0.5 Jim in the largest size channels. The instrument is 
calibrated prior to use using monodispersed spheres. 

( Passive Cavity Aerosol Spectrometer Probe (Airborne) PMS 
Model PCASP-100X 0.10 - 3.0 mm Operating Manual) 
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5. Particle Soot Absorption Photometer (PSAP) 

The PSAP measures the aerosol absorption coefficient of 
an air sample based on the integrating plate technique. An 
external vacuum source provides an air sample to the PSAP; 
differences in the optical transmission of a filter due to 
aerosol deposition are converted to absorption using Beer's 
Law and a calibration transfer coefficient determined by the 
filter type. The PSAP is designed to give a continuous 
measurement of absorption coefficient with time averaged 
resolution of 5 to 300 seconds. The absorption coefficient 
is measured at 0.565 mm with a sensitivity of lO'Si' 1 for 1 
minute averages. ( Particle Soot /Absorption Photometer (PSAP) 
Operation Procedures ) 
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IV 



OPTICAL DEPTH RETREIVAL PROCEDURES 



This chapter outlines the procedures used to calculate 
aerosol optical depths from the satellite imagery included 
in the reference data sets. 

A. SATELLITE IMAGE RETREIVAL/ DISPLAY 

Satellite image retrieval and display was performed on 
the Terascan Earth Remote Sensing System by SeaSpace 
Corporation. Terascan allows display and enhancement of 
NOAA POES High Resolution Picture Transmission (HRPT) data 
and GOES variable (GVAR) formatted data at full resolution. 
Embedded routines in the Terascan system convert sensor 
radiance counts into albedo/brightness temperature, 
calculate required angles from Sun-Earth-satellite geometry, 
and map imagery from different sources to the same spatial 
grid for comparison. Imagery used in this study was 
previewed in Terascan prior to optical depth retrieval to 
conduct cloud and Sun glint screening. Terascan also offers 
a full suite of post processing enhancements which greatly 
aided the optical depth analysis. 

B. OPTICAL DEPTH RETRIEVAL 

The optical depth retrieval technique used for the GOES 
Imager and NOAA AVHRR data is an automated process which 
required the use of both Terascan and FORTRAN 77 code. The 
processing sequence is illustrated Figure 4.1. 

Pre-processing consists of the use of Terascan commands 
to retrieve raw satellite images from the pass disk of the 
Terascan receiver system. Manual cloud and Sun glint 
screening is performed in this step. 

The main processing of the data was completed using a 
mixture of Unix Borne scripts (with embedded Terascan 
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commands) and compiled FORTRAN code. The Unix scripts 
registered the raw files to a predetermined geographic area, 
calibrated the data (if required), calculated significant 
angles (satellite zenith. Sun zenith, relative azimuth, 
and scattering angle) , and created an export binary file for 
mathematical processing in FORTRAN. The FORTRAN code 
calculated all radiative transfer calculations including 
optical depth. After calculation, the radiative transfer 
output from the FORTRAN code was imported back into Terascan 
format and assembled with the initial channel data. 

Post-processing included the use of Terascan display 
and enhancement features to analyze the imagery. 

C. RADIATIVE TRANSFER CODE 

The FORTRAN code used to calculate optical depth from 
measured satellite radiance follows the theory of Chapter II 
and calculates a solution to Equation 4.10. The central 
value of the satellite channel's spectral response is used 
as the input wavelength. Calculation of the specific input 
variables into Equation 2.6 are described in this section. 

1. Solar Irradiance/Solar Radiance 

Input values for solar irradiance are determined by 
calculating a weighted average of solar irradiance 
integrated across the satellite radiometer channel spectral 
response. For the NOAA 12 and 14, values were taken from 
the NOAA Polar Orbiter Data Users Guide (Kidwell, 1995) . 
For the GOES, spectral response functions for the GOES 8 and 
GOES 9 documented in the GOES-IJ/KLM SN03 Imager Data and 
Calibration Handbook (1994) and GOES Calibration and 
Alignment Handbook for the Imager SN04 Instrument (1994) 
were used to calculate the weighted solar irradiance. Table 
4.1 lists the values of solar irradiance used for the 
various satellites/channels. Solar radiance is calculated 
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by dividing solar irradiance by pi (E 0 /7i) . Prior to input 
into the model, E 0 is corrected for variation in the Earth- 
Sun distance. 



Satellite 


Radiometer Channel 


Eo (W/m 2 pm) 


NOAA 12 


AVHRR Ch 1 


1614 


NOAA 12 


AVHRR Ch 2 


1050 


NOAA 14 


AVHRR Ch 1 


1628 


NOAA 14 


AVHRR Ch 2 


1030 


GOES 8 


Imager Ch 1 


1629 


GOES 9 


Imager Ch 1 


1617 



Table 4.1. Values of Solar Irradiance 



2. Ozone and Rayleigh Optical Depths 

Input values for ozone and Rayleigh optical depths are 
determined by calculating a weighted average of reported 
ozone and Rayleigh optical depths by Elterman (1970) 
integrated across the satellite radiometer channel spectral 
response. Table 4.2 lists the values of both ozone and 
Rayleigh optical depths used for the various 
satellites/channels . 



Satellite 


Radiometer Channel 


Ozone 8 


Rayleigh 8 


NOAA 12/14 


AVHRR Ch 1 


0.027 


0.057 


NOAA 12/14 


AVHRR Ch 2 


0.0021 


0.019 


. GOES 8/9 


Imager Ch 1 


0.024 


0.061 



Table 4.2. Values of Ozone and Rayleigh Optical Depths 



3. Scattering Phase Function 

Parameterization of the scattering phase function is 
the most difficult part of the optical retrieval. Knowledge 
of aerosol size distribution is required to accurately 
calculate the phase function from Mie theory. Since this 
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distribution is not known exactly without in situ aerosol 
measurement, the scattering phase function must be 
parameterized by other methods. For this study, two methods 
are used. One method uses an empirically derived aerosol 
size distribution to calculate the scattering phase 
function. The second method uses multispectral radiance 
differences measured by a satellite to parameterize the 
scattering phase function. These methods are described 
below. 



a. NOAA/NESDIS derived Scattering Phase Function 
An empirically derived aerosol size distribution 
provided by the NOAA/NESDIS Satellite Research Laboratory is 
used to calculate the scattering phase function. The 
aerosol size distribution is a single-mode log-normal 

distribution with a mode radius of 0.1 |im, standard 

deviation (a) of 2.03, and a complex refractive index m = 
1.4+0.0i (Ignatov et al, 1995). This type of scattering 
phase function is applied statically across the entire image 
during retrieval and does not account for variations in the 
aerosol size distributions known to exist, especially in 
coastal regions. Figure 4.2 and 4.3 illustrate the aerosol 

size distribution and phase functions at 0.63 |im and 0.86 

Hm. 



b. Particle Size Parameter (S 13 ) Derived 
Scattering Phase Functions 

Durkee et al . (1991) proposed a method of 

parameterizing the scattering phase function based on a 
ratio of the aerosol radiance measured in the channels 1 
(visible) and 2 (NIR) of the NOAA AVHRR. Because scattering 
efficiency (Q scaC ) of an aerosol distribution is wavelength 
dependent, scattering for a specific aerosol population 
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peaks when the radius of the aerosol is nearly equal to the 
radiation wavelength. Subsequently, radiance counts 
measured by the AVHRR visible and NIR channels will change 
with aerosol size distribution changes such that the ratio 
of channel radiances will be larger for smaller size 
particle distributions and smaller for larger size particle 
distributions . 

Durkee et al . (1991) termed the ratio of the 
channel aerosol radiances the particle size parameter, S 12 . 
Since S 12 varies pixel-by-pixel across the entire satellite 
image, the scattering phase functions can be parameterized 
pixel-by-pixel, allowing variations in aerosol distributions 
to be properly factored into the optical depth retrieval. 

For the optical depth retrieval, the scattering 
phase functions and extinctions for seven model aerosol size 
distributions (M0-M6) are calculated using Mie theory. 
These distributions consist of one single-mode and six two- 
mode log-normal distributions with varying mode radii and 
standard deviations designed to model the typical variations 
of aerosol distributions in the marine atmosphere. The 
first mode models the background aerosol while the second 
mode models ocean-produced aerosol. The NOAA/NESDIS aerosol 
size distribution represents a rough approximation to the 
average of these distributions. The refractive index used 
is the same as the NOAA/NESDIS model. Table 4.3 lists the 
mode radii and standard deviations used. Figures 4.4, 4.5, 
and 4.6 illustrate the resulting aerosol size distributions 

and the corresponding scattering phase functions at 0.63 |im 
and 0.86 nm. 
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Model 


Mode Radii (nm) 


Std Dev (a) 


MO 


0. 1/0.0 


1. 7/0.0 


Ml 


0.1/0. 3 


1.7/2. 1 


M2 


0.1/0. 3 


1.7/2. 2 


M3 


0.1/0. 3 


1.7/2.35 


M4 


0. 1/0.3 


1.7/2.51 


M5 


0.1/0. 3 


1.7/2. 6 


M6 


0. 1/0.3 


1. 7/2.7 



Table 4.3. Mode Radii and Standard 
Deviations for Model Aerosol Size 
Distributions 



The calculated scattering phase functions and 
extinction coefficients' for these distributions are 
converted to S 12 using the following equation (based on 
linear single scattering theory) : 



^12 




Pchl^ochl* 7 extchl 
Pch2 F o c h2 a extch2 



Due to water vapor absorption in channel 2 of the AVHRR, a 
correction to S 12 is made following Mahony (1991) based on 
the split-channel (channels 4 and 5) water vapor retrieval 
proposed by Dalu (1986) . The resulting S 12 values for each 
aerosol model are shown in Figure 4.7. 

During processing of the satellite radiance data 
in the optical depth retrieval code, S 12 values are 
calculated for each pixel from AVHRR channel 1 and channel 2 
data. Computed S 12 values along with scattering angle are 
then entered into a lookup table representing Figure 4.7 to 
determine the model aerosol distribution (Figure 4.4) best 
represented by the observed radiance measurements. Using 
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the scattering angle and model size distribution, phase 
function values are then selected from lookup tables based 
on Figures 4.5 and 4.6. 
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Figure 4.1. Satellite Aerosol Optical Depth Retrieval. 
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NOAA Model Aerosol Size Distribution 




Figure 4.2. NOAA/NESDIS Aerosol Size Distribution. 



NOAA Model Phase Functions 




Figure 4.3. NOAA/NESDIS Scattering Phase Functions. 
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Model Aerosol Size Distributions 




Figure 4.4. Model Aerosol Size Distributions. M0-M6 
correspond to model two-mode, lognormal aerosol 
distributions with MO representing the background 
(continental) aerosol mode and M6 representing the largest 
oceanic aerosol mode. 



Model Phase Functions at .63 nm 




Figure 4.5. Model Phase Functions at 0.63 |jm. M0-M6 
correspond to model two-mode, lognormal aerosol 
distributions with MO representing the background 
(continental) aerosol mode and M6 representing the largest 
oceanic aerosol mode. 
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Model Phase Functions at .86 pm 




Figure 4.6. Model Phase Functions at 0.86 |lm. M0-M6 
correspond to model two-mode, lognormal aerosol 
distributions with MO representing the background 
(continental) aerosol mode and M6 representing the largest 
oceanic aerosol mode. 



Model Sj 2 Values 




Figure 4.7. Model S 12 Values. M0-M6 correspond to model two- 
mode, lognormal aerosol distributions with MO representing 
the background (continental) aerosol mode and M6 
representing the largest oceanic aerosol mode. S 12 values 
are calculated using linearized, single-scattering theory. 
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V 



RESULTS 



A. OVERVIEW OF TARFOX 

During the TARFOX period, 10-31 July 1996, coordinated 
aircraft measurements and satellite overpasses were limited 
due to a number of factors. The weather played a major role 
due to a series of low pressure systems and associated 
cloudiness repeating approximately every 2-3 days over the 
TARFOX area. Another limitation was the viewing geometry of 
the NOAA 14, eliminating 6 NOAA 14 overpasses from the data 
set due to Sun glint contamination (10, 17, 20, 27, 28, 29 
July) . Based on all factors, the data set selected for this 
thesis included aircraft, NOAA 14, and GOES 8 data on 16, 
23, and 25 July. Of these days, only 25 July had a complete 
set of spatially and temporally coordinated sunphotometer , 
PCASP, PSAP and satellite data (partially due to data access 
limitations imposed by some of the TARFOX participants) . 
Additionally, it was the most cloud-free day analyzed (see 
Figure 5.1). Subsequently, the majority of the results of 
this study are concentrated on data from 25 July. An 
additional day of coordinated sunphotometer data from a land 
site on Bermuda and NOAA 14 overpass on 18 July are also 
analyzed; unfortunately, GOES 8 imagery extending over 
Bermuda was not archived during this period due to storage 
limitations. Some UK C-130 PCASP data taken on days 
without useful NOAA 14 imagery are also presented to 
illustrate variability of aerosol loading in the TARFOX 
area. Table 5.1 lists the aircraft flights, satellite 
overpasses and data used. 
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Date 


Time 


Aircraft/ Satellite 


Data 


16 July 96 


1840 UTC 


NOAA 14 


AOD 


16 July 96 


1853 UTC 


UW C-131A 
(Flight 1727) 


sunphotometer 


18 July 96 


1816 UTC 


NOAA 14 


AOD 


18 July 96 


1656- 
1956 UTC 


Bermuda 


sunphotometer 


21 July 96 


1414 UTC 


UK C-130 
(Flight A466) 


PC ASP 
PSAP 


23 July 96 


1906 UTC 


NOAA 14 


AOD 


23 July 96 


1839 UTC 


UW C-131A 
(Flight 1732) 


sunphotometer 


25 July 96 


1432 UTC 


GOES 8 


AOD 


25 July 96 


1435 UTC 


UW C-131A 
(Flight 1734) 


PCASP 


25 July 96 


1842 UTC 
1906 UTC 


UW C-131A 
(Flight 1735) 


sunphotometer 

PCASP 


25 July 96 


1842 UTC 


NOAA 14 


AOD 


25 July 96 


1845 UTC 


GOES 8 


AOD 


25 July 96 


2050 UTC 


UK C-130 
(Flight A469) 


PCASP 

PSAP 


25 July 96 


2045 UTC 


GOES 8 


AOD 


27 July 96 


2203 UTC 


UK C-130 
(Flight A470) 


PCASP 

PSAP 



Table 5.1. TARFOX Data and Sensors Used in Study. 



B. 25 JULY 1996 
1. Weather 

Pre-frontal conditions existed over the TARFOX area on 
25 July 1996 with southeasterly flow at 5-10 kts due to a 
1015 mb low pressure center forming inland over the North 
Carolina/Virginia border. Figure 5.1 provides a NOAA 14 
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AVHRR channel 1 visible image of the experiment area. A 
broad area of low clouds associated with an easterly moving 
low pressure center is present in the northeast corner of 
Figure 5.1. Some cumuloform cloudiness associated this low 
center's dissipating trailing cold front exists in the 
southern region of the image. Otherwise, generally clear 
conditions were present in the flight areas near 38N 74. 5W. 
With the southeasterly flow, an oceanic aerosol distribution 
in expected off the coast of the U.S.. 

2 . Absorption 

In order to test the assumption that absorption due to 
aerosols in the cloud-free marine environment is much 
smaller than scattering and, therefore, negligible, PSAP 
data measured by the UK C-130 was analyzed. Figure 5.2 
provides the PSAP measured absorption coefficient for flight 
A469 on 25 July 1996. Additionally, PSAP absorption 
coefficients for flights A466 and A470 are included for 
reference. Comparison of Mie calculated scattering 
coefficients based on PCASP aerosol size distributions for 
the same flight profiles suggest that absorption is at least 
an order of magnitude smaller than scattering. Resulting 

estimates show that co 0 * 0.97-0.99. Therefore, assuming 

zero absorption results errors in co 0 of less than 5%. At 

0.63 pm, co 0 = 0.98 corresponds to a complex index of 
refraction of 1.4 + 0.002i. Using this value in phase 
function calculations based on PCASP size distributions from 
flight A469 indicate that assuming zero absorption also 
leads to less than 5% error in the MIE calculated phase 
function (see Figure 5.3). 
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3 . Water Vapor 

In order to validate water vapor corrections in channel 
2 of the AVHRR, estimates of column water vapor were made 
using temperature and dew point measurements during flight 
1735. Figure 5.4 provides the calculated water vapor 
profile measured during an ascent of flight 1735. In order 
to estimate total column water vapor, it was assumed that 
water vapor above the maximum height of the profile was zero 
and a 4 th order polynomial fit to the profile was integrated 
across the depth of the profile. This resulted in a value 
of 1.77 gm/cm 2 . Using the Dalu (1986) split-channel (AVHRR 
channel 4 - channel 5) technique in a 10x10 pixel box over 
the ascent region of flight 1735, total column water vapor 
was retrieved from the 1842 UTC NOAA 14 overpass. 
Resulting values of retrieved total column water vapor from 
the AVHRR were 2.59 gm/cm2 with a standard deviation of 
0.14. Both values lie within the climatological extremes of 
atmospheric water vapor; uncertainty in the amount of water 
vapor above the measured profile possibly accounts for the 
differences in values. Comparison of the effects of these 
differences on S 12 and retrieved phase function indicate 
that a 3 0% difference in water vapor results in an 
approximate difference of 1-2% in S 12 and less than 3% in 
phase function. 

4. Aerosol Distributions/Phase Functions 
a. 25 July Intercomparisons 

Comparison of PSCAP data collected during 
ascent/decent of flights 1734, 1735, and 1736 to model 
aerosol size distribution provided some interesting results. 
Note that absolute calibration between the PCASP instruments 
used on the two aircraft had not been accomplished at the 
time of this study; therefore, some differences between 
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measured aerosol distributions may be due to the use of two 
instruments. Figures 5.5 and 5.6 contain measured and model 
aerosol size distributions in terms of number concentration 
(dN/dr) and volume (dV/dlogr) . Good correlation between the 
aircraft measured distributions is illustrated, suggesting a 
nearly homogeneous air mass present over the TAFROX area 
both spatially and temporally. Based on the measured S 12 
values from the 1842 UTC NOAA 14 overpass, phase functions 
based on aerosol models M2 and M3 were selected by the 
satellite retrieval algorithm and are included for 
comparison. 

Significant departures in both dN/dr and dV/dlogr 
between measured and model distributions in the 0.1 to 0.7 

pm radius range are noted. At radii larger than 0.7 pm, 
models M2 and M3 fit the measured data well in dN/dr. 
Distribution differences are more apparent in the dV/dlogr 
plots; however, the two-mode M2 and M3 models demonstrate 
some skill at modeling the mode radii of the measured 
distributions but not the amplitude (up to 95% error) . A 

third, smaller mode near 0.8 pm is not modeled in M2 or M3. 
The NOAA/NESDIS single-mode model is the worst fit to the 
measured distributions; however, this distribution broadly 

captures the mode near 0.8 pm in both dN/dr and dV/dlogr. 

Figures 5.7 and 5.8 profile the potential 
temperature and relative humidity (RH) observed along the 
ascent/descent tracks of flights 1735 and A469. Note the 
consistency in the lower tropospheric structure. The 
shallow mixed layer (250m) and relatively high RH (50-60%) 
between 250 m and 3500 m is indicative of a maritime air 
mass. This is consistent with the wind field observed on 25 
July 1996. The characteristically high number and low 

volume of particles < 0.7 pm observed may be due to the 
longevity of the aerosol suspended over the ocean. 
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Figures 5.9 and 5.10 show the calculated phase 
functions from both measured and modeled size distributions 

at 0.63 |im and 0.86 pm, respectively. With the exception of 

the NOAA/NESDIS phase function at 0.86 |im, all modeled phase 
functions are 20-40% larger than the PCASP derived phase 
functions in the 160-180° back scattering angles (scattering 

angle for 1842 UTC NOAA 14 overpass 168°) . Despite errors 
in magnitude, M2 and M3 fit the shape and slope of the 
derived phase functions. Characteristic sinusoidal patterns 
in the NOAA/NESDIS phase function are not present in the 
measured distribution phase functions. Also, in comparison 

of PCASP derived phase functions at 0.63 Jim and 0.8 6 Jim, an 
increase in phase function across all scattering angles is 
observed between the two wavelengths. This characteristic 
is not observed in the NOAA/NESDIS phase function, but is 
well represented in the M1-M6 phase functions. 

Based on these observed differences in measured 
and modeled phase functions on 25 July 1996, errors in 
satellite retrieved aerosol optical depth should be on the 
same order due to the linearity of the algorithm (Equation 
2 . 6 ) . 



b. 21, 25, and 27 July Intercomparisons 
To understand the variability of the aerosol 
distributions during TARFOX, variations in UK C-130 PCASP 
measured aerosol size distributions, atmospheric boundary 
layer profiles, and calculated phase functions collected 
during flights A466, A469, and A470 are analyzed. Figures 

5.11 and 5.12 provide comparisons of observed and modeled 
distributions in number concentration and volume, 

respectively. Figures 5.13 and 5.14, the corresponding 

boundary layer profiles, show a definite change in air mass 
occurs between 21 July and 25 July and again from 25 July to 
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27 July. The depth of the mixed layer (~ 1500 m) and the 
low relative humidities above this level on 21 and 27 July 
suggest an air mass of continental origin. This is 
consistent with northwesterly flow observed over the TARFOX 
area on 21 and 27 July. With the change in air mass, large 
variations in both aerosol number concentration and volume 
are present, especially in the smallest mode. Variability 
in the aerosol distribution is expected; however, modeled 
variability does not adequately capture the observed aerosol 
fluctuations. Because of lower aerosol number concentration 
and volume on 21 and 27 July, model estimates between MO and 
M3 are expected over the entire period. 

From Figures 5.15 and 5.16, error in the model 
size distributions lead to approximately 20-40% error in 
the model phase function based on expected S 12 . These 
results are consistent with the results of the 
intercomparison of data on 25 July. 

5. NOAA 14 Aerosol Optical Depth Retrieval 

Sunphotometer aerosol optical depth data collected 
onboard flight 1735 during low-level portion of the flight 
was used in the validation of NOAA 14 retrieved aerosol 
optical depth from the 1842 UTC overpass. Figure 5.17 
demonstrates good agreement between sunphotometer measured 
extinction and extinction calculated from measured size 
distributions from flights 1734, 1735 and A469, especially 
at NOAA 14 channel 1 and 2 wavelengths . The sunphotometer 
data represented in Figure 5.17 reflects an average of the 

low level portion of flight 1735 with error of ± 6% based on 
standard deviation of the measurements. Figure 5.18 
presents sunphotometer aerosol optical depths collected 

during flight 1735. Only data at 0.38 (im, 0.451 urn, 0.525 
pm and 1.021 pm were recorded during the flight. Values at 
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0.63 pm and 0.86 pm were derived from a power law fit to the 
reported data. 

Figure 5.19 represents the 1842 UTC NOAA 14 retrieved 
aerosol optical depth; the low level portion of flight 1735 
is represented by the red track. Figures 5.20 and 5.21 
represent comparisons of satellite retrieved aerosol optical 
depth collected in a box surrounding the flight 1735 flight 

track and sunphotometer data at 0.63 pm and 0.86 pm, 
respectively. Despite large errors in the phase function 
reported above, the satellite data fits closely to the 
sunphotometer observed data in magnitude and shape. This 
result suggests that the retrieval technique is sensitive to 
observed aerosol variation within the sample region. 

A multiple-scattering radiative transfer model with bi- 
directional surface reflectance was used to investigate this 
unexpected fit. The model, the Second Simulation of 
Satellite Signal in the Solar Spectrum (Six S) model 
(Vermote et al . 1995), includes surface reflectance due to 
direct and diffuse radiance. After conducting simulations 
using observed PCASP aerosol distributions from flight 1735, 
observed Sun-Earth-satellite geometries along the flight 
1735 track, and observed wind speed/direction, analysis of 
sunphotometer. Six S, and NOAA 14 AVHRR aerosol optical 
depths were completed. The calculated Six S model aerosol 

optical depth at 0.63 pm fit the sunphotometer observed data 
to 0.01 optical depth accuracy. Comparison of Six S 
calculated aerosol radiance (8.815 W/m 2 sr pm) to NOAA 14 

derived aerosol radiance (9.11 W/m 2 sr pm) showed close 
agreement; this suggests that linear treatment of Rayleigh 
scatter and foam and subsurface reflectance in the linear 
single scatter model is comparable to the non-linear Six S 
model. Figure 5.22 graphically illustrates the overall 
results of the comparison of the models. Four possible 
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solutions to the problem are represented using observed 
aerosol distributions: (1) Six S with a full bi-directional 
surface reflectance model, (2) Six S without a bi- 
directional surface reflectance model (multiple scattering 
solution), (3) linear single scattering solution without a 
bi-directional surface reflectance model, and (4) the linear 
single scattering solution with the bi-directional surface 
reflectance approximation used in this study. In a linear 
sense, the slope of these lines is representative of the 

effective phase function, where L a = PefftVs)^* Error 
between the multiple scattering solution without a specular 
surface and the linear single scattering solution without a 
specular surface is approximately 17% for aerosol optical 
depths observed. Error between the full Six S solution and 
the linear single scattering solution used in this study is 
approximately 22%. Since the Six S solution is non-linear', 
it is not possible to separate the contribution to satellite 
observed radiance due to multiple scattering from that due 
to spectral surface reflectance. Therefore, the total error 
in the satellite retrieved aerosol optical depths using the 
study algorithm is due in part to both lack of multiple 
scattering and approximating spectral surface reflection. 
By chance, the error in the phase function representation in 
the model is offset by these errors, resulting in the fit of 
retrieved aerosol optical depth to the sunphotometer data. 

6. GOES 8 Aerosol Optical Depth Retrievals 

Despite the compensating errors in the satellite 
aerosol optical depth retrieval algorithm, the algorithm was 
extended to the 1845 UTC GOES overpass. Figure 5.23 
provides the retrieved GOES 8 aerosol optical depths using 
the 184 2 UTC NOAA 14 S 12 phase function parameterization. 
Figure 5.24 demonstrates the fit of the GOES data to the 
sunphotometer data, again in good agreement. 
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Based on the consistency of the aerosol distributions 
taken during flights 1734, 1735, and A469, S 12 data from the 
1842 UTC NOAA 14 overpass were applied to GOES 8 imagery at 
1432 UTC and 2045 UTC. Figures 5.25 and 5.26 provide 
histograms of calculated aerosol optical depth and phase 
functions for these passes as well as the 1842 UTC NOAA 14 
and 184 5 UTC GOES 8 overpasses for an area of the imagery 
which remained cloud-free over the time span of the passes. 

The similarity of the GOES 8 and NOAA 14 aerosol 
optical depth data in Figure 5.25 demonstrates the potential 
of the GOES Imager data in aerosol optical depth retrieval. 
Throughout the time period from 1432 to 2045 on 25 July 

1996, the GOES 8 Imager scattering angle varied from 120° to 

145° while the NOAA 14 ■ AVHRR scattering angle was 168°. 
Despite the variation in scattering angle and subsequent 
variation in phase function, consistency in retrieved 
aerosol optical depths between the two instruments was 
observed (Figure 5.25). Although the error in magnitude of 
the model phase function was on the order of 2 0-4 0%, the 
overall shape and slope of modeled phase functions had to be 
well represented in the algorithm in order to account for 
the close fit of the retrieved aerosol optical depths. Due 
to the linearity of the solution, any significant deviation 
in the shape and/or slope of the model aerosol distribution 
to the observed aerosol distribution would have resulted in 
significant variation in retrieved optical depths at each 
scattering angle. In the homogeneous offshore environment 
on 25 July, application of S 12 forward and backward in time 
appears reasonable. 

It is expected that degradation of the S 12 information 
with time will occur and S 12 application to GOES data will 
be limited temporally. These limitations require further 
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study with a larger data set in varied aerosol loading 
conditions. 

Comparison of Figures 5.20 and 5.24 also illustrate the 
aerosol optical depth resolution expected from the two 
instruments based on digital resolution; the banding in the 
satellite data is due to the digitization steps of the 
satellite sensor. Aerosol optical depth resolution of the 
NOAA 14 AVHRR appears to be on the order of 0.04, while GOES 
8 resolution is on the order of 0.8. 

C. 18 JULY 1996 

1. Weather 

The focus area for data on 18 July 1996 was near 
Bermuda. The synoptic pattern over the island was 
characterized by a 1026 mb high resulting in generally clear 
skies with small clusters of low level cumulus. Winds 5-10 
kts from the southwest were observed over Bermuda. Due to 
the trajectory of the air mass from the southwest, aerosol 
loading over Bermuda was relatively low due to lack of 
continental influence. 

2. NOAA 14 Aerosol Optical Depth Retrieval 

This case was selected for analysis based on low 
optical depths observed in the land based sunphotometer data 
over Bermuda. Figure 5.27 provides the histogram comparison 

of sunphotometer aerosol optical depth data (at .699 |im) 
taken between 1659 and 1959 UTC 18 July 1996 on Bermuda to 
aerosol optical depth retrieved from 1816 UTC NOAA 14 
channel 1 in a cloud-free area just to the southwest of 
Bermuda. Note the good agreement in the data. The 
algorithm, despite noted errors, correctly computed lower 
aerosol optical depths. Also, this example shows the 
expected resolution of the NOAA 14 retrieved aerosol optical 
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